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Abstract. Infrared observations of hot massive stars and their environments provide a detailed 
picture of mass loss histories, dust formation, and dynamical interactions with the local stellar 
medium that can be unique to the thermal regime. We have acquired new infrared spectroscopy 
and imaging with the sensitive instruments onboard the Spitzer Space Telescope in guaranteed 
and open time programs comprised of some of the best known examples of hot stars with cir- 
cumstellar nebulae, supplementing with unpublished Infrared Space Observatory spectroscopy. 
Here we present highlights of our work on the environment around the extreme P Cygni-type 
star HDE316285, providing some defining characteristics of the star's evolution and interactions 
with the ISM at unprecented detail in the infrared. 
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1. Introduction 

Observations of the circumstellar environments of Lumnious Blue Variables (LBVs) 
with ESA's Infrared Space Observatory (ISO) and modern ground-based imaging devices 
have demonstrated how sensitive mid-infrared imaging and spectroscopy can substan- 
tially improve our knowledge of the distribution of material, the physical and chemical 
properties of nebular dust and gas, and therefore on the history of mass loss from the 
surface of the massive central star during its evolution from the Main Sequence. Other 
than the bright and compact environments of most LBV nebulae, none of the lower 
surface brightness nebulae around OB or Wolf-Rayet stars could be observed with ISO 
due mainly to instrumental sensitivity limitations. The instruments on the Spitzer Space 
Telescope have since offered significant gains in sensitivity and sky coverage; for com- 
parison, the Spitzer Infrared Spectrometer (IRS; 5.3 - 38 /im) is factor of ~50-100 more 
sensitive than the ISO Short Wavelength Spectrograph (SWS; 2.4-45.2 /im) at 5 /im, 
depending on resolution modes of the two instruments. We have exploited the Spitzer 
IRS, the Infrared Array Camera (IRAC) with imaging bands at 3.4, 4.6, 5.8, and 8.0 /im, 
and the Multiband Imaging Photometer for Spitzer (MIPS) with imaging capabilities at 
24, 70, and 160 /im in Guaranteed and Open Time programs to observe the environ- 
ments around a number of hot, massive stars. Most of these stars are surrounded by ring 
nebulae and are well known for their optical and/or radio properties; NGC2359 (WR), 
Ml-67 (WR), G79.29+00.46 (Ble), and HD148937 (Ofp?e) are examples. 

In this paper we summarize our study of the massive B supergiant and candidate LBV 
HDE316285 from our program. This highly luminous P Cyg-type star has been suspected 
of being surrounded by an extended, cold nebula (McGregor, Hyland, & Hillier 1998); 
we have obtained conclusive observations of the environment around this star, which 
harbor telling characteristics of the star's mass loss history, wind-wind and wind-ISM 
interactions, the molecular content of the gas and thermal history of the dust. 
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2. The ring nebula around the extreme P Cygni star HDE316285 

The combined 1.4-38 /am spectrum of HDE316285 is shown in Fig. [H This star has 
been quantitatively characterized by Hiher et al. (1998) using hne-blanketed non-LTE 
wind models as an "extreme" P Cygni-type star, with an optical and near-IR stellar wind 
spectrum that is quite similar to that of P Cyg. Spectral variability, the stellar properties 
and chemical content of HDE 316285 are similar to known LBVs, but is more extreme 
than P Cyg because of its high wind performance number (= ratio of wind momentum 
to radiative momentum), some 30 times greater. This number, however, is based on the 
assumption of a smooth, homogeneous wind and is subject to downward revision in the 
case of a clumped wind and associated reduction to the mass loss rate (cf. Hillier&Miller 
1999). Somewhat unexpectedly, however, the infrared continuum, which can be used to 
constrain the volume filling factor (e.g., Morris et al. 2000), is contaminated by a strong 
thermal excess which becomes evident at around 16/im (see Fig. [T]). In fact McGregor, 
Hyland, & Hillier (1988) had pointed out that HDE316285 is a moderate IRAS source, 
exhibiting cold dust (~60 K) removed from the central star. Since this star is projected 
^ 1° from the Galactic Center and is now recognized to be in close (apparent) proximity 
to the infrared star forming region (SFR) Sgr D (discussed below), a significant thermal 
component arising from surround molecular material heated by nearby young stars would 
provide the most likely explanation for the excess emission in the IRAS 25 and 60 /im 
passbands. 

Infrared imaging. We now confirm that IIDE316285 is definitely surrounded by its own 
circumstellar nebula. Following careful inspection of archival Midcourse Space Experi- 
ment 21.3 /im data in which a bipolar nebula can be seen just barely above the detection 
threshhold, our MIPS 24/im imaging (see Fig. [2]) reveals a limb-brightened, 3'.4x4M ring 
nebula with the major axis oriented ~ 40° from the north. The nebula is not discernable 
overall in the IRAC bands, though one bright knot of 8.0 /im emission is detected ~ V 
NE of the central star. This knot and the brightened southern limb (respectively labelled 
"2" and "1" in Fig. [l]) exhibit peak nebular fluxes of 624 and 450 MJy/sr, and were 
targeted by us for IRS spectroscopy, shown in Fig. [2] and discussed below. The nebula at 
24 /im also clearly exhibits a clumpy structure, and filaments directed radially outwards 
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Figure 1. Combined near-IR (1.4-2.4/im), ISO/SWS (2.4-10/im), and Spitzer/IRS (10.0-38/im) 
SEDs of HDE316285 (upper) and G79.29+00.46 (lower). From Morris et al. (2008) submitted. 
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Figure 2. Left: HDE316285 at 24 /xm with Spitzer/MIPS. The image is 5'. 5 x 8, north is up, 
east is left. Two regions targeted for spectroscopy are indicated. Right: IRS 5.3-38/im spectra 
at two positions in the nebula (background-subtracted), with main features identified. 

particularly evident in the NE half. The SW rim of the nebula appears to be interacting 
with material emitting at 8.0 fim in the surrounding ISM. (see also Fig. [3]). 

Spectral content: PAHs and rotational H2. Spectra extracted at the two positions indi- 
cated in Fig. [2] and fully corrected for background emission with a dedicated off-position 
observation reveal thermal spectra with two principle components: a dominant cool com- 
ponent in the 75-125 K range, and a warm component at Tdust > 350 K. Atop the warm 
component we see the family of polyaromatic hydrocarbon (PAH) bands at 6.22, 7.63, 
8.63, 11.22, and 12.75 /im, and possibly a number of secondary peaks. These bands are 
common in reflection nebulae and the shells of AGB and post-AGB stars with mixed C 
and O chemistries, but they have also been detected mixed with the crystalline silicate 
dust in the nebulae around LBVs AG Car and R71 (Voors et al. 1999, 2000). We return 
to these features at the end of this paper. 

The HDE316285 nebular spectra also exhibit pure rotational lines of H2 arising from 
optically thin quadrupole transitions, from S(0) at 28.3 /im through S(6) at 6.1 /im, 
though the S(4) 8.0 jam and S(6) lines are blended with PAH emission. This is the first 
known detection of rotational transitions of H2 in the nebula of a hot massive star. St. 
Louis et al. (1998) have previously detected the ro-vibrational 1-0 S(l) 2.112 /im line 
in the NGC2359 nebula around WR7, in a region where the nebula may be interacting 
with a surrounding bow shock. From this single line, however, St. Louis et al. could 
not deduce the excitation mechanism. This is relavent to the debate on the contribution 
of hot stars to the H2 luminosity in starburst galaxies and ULIRGs, where flourescent 
excitation in the UV radiation fields of OB and WR stars has been generally preferred 
over collisionally-induced emission for lack of cases of the latter in local settings. 

The H2 rotational lines are readily thermalized at moderate volume densities, allowing 
us to estimate H2 column densities and excitation temperatures under the assumption 
of LTE in order to gain insight into the excitation conditions of the line forming regions. 
We follow excitation diagram methods developed by Burton (1992) and Gredel (1994), 
and applied to comparable environments (as in, e.g., NGC7129 by Morris et al. 2004 and 
references therein). The column density for transition j is Nj = AttFj /QEjgjAj, where Ej 
is the energy of the upper level, gj is the statistical weight of ortho- and para-transitions, 
and Aj is the transition probability. The line column density is related to the total column 
density and excitation temperature Tex by \n{Nj/gj) = ln{Ntot/Q) — Ej / {kT^^)^ and 
may be solved for with a Boltzmann excitation diagram. The partition function Q was 
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Figure 3. Mosaic of the Sgr D region, with radio sources identified by Mehringer et al. (1998) 
indicated. Blue corresponds to 4.6 /xm, green to 8.0 /xm, and red to 24 /xm. 

determined with an ortho-para ratio of 3 under the assumption of thermal equilibrium. 
For the observed line fluxes per unit solid angle Fj/Q^ we we have extracted over regions 
of 40 arcsec^ and 62 arcsec^ in the 5.3 - 14 /im and 14 - 35 /im ranges, respectively. In 
these areas we must view our measurements as spatially integrated. The measured line 
intensities were corrected for reddening Ay = 6.0 mag (Hillier et al. 1998). Measurement 
uncertainties of the S(3) - S(5) lines arise predominantly from the blending with the PAH 
features and the flux calibration. The S(6)/S(5) or S(4)/S(5) line ratios could be also 
affected by an ortho-para ratio different than the canonical (thermal equilibrium) value 
of 3, as in the presence of shocks, or by the superposition of multiple layers along the 
line of sight with different physical conditions. From our computed values of \n{Nj/gj) 
and E^p/K we obtain estimates of the total H2 column density A/tot — 1.5 x lO^^cm"^ 
and excitation temperature Tex — 240 K, averaged over the two positions. 

The relatively low value of Tex at the high H2 column density (compared to typical 
photodissociation regions) seem quite reasonable for line emission in the outflow arising 
from collisional (dissociative J-shock) excitation. We favor the J-shock mechanism over 
C-shocks in the outflow, based on the S(0) line intensity and strong [Si ll] 34.815 jam 
emission. From J-shock models (Hollanbach & McKee 1989) we see also by the presence 
of [Ne 11] 12.81 /im that the velocity is probably greater than 50 km sec~^. 

Is HDE316285 a hyper-luminous counterpart of the Pistol Star? In order to estimate 
the physical size and age of the nebula around HDE316285 we must know the distance 
to the star. Hillier et al. (1998) adopted 2 kpc, as a rough average of literature values 
between 1 and 3.4 kpc, and conceded that there is no real constraint on the distance. 
At that time it was not recognized that HDE316285 lies in close (apparent) proximity to 
the infrared SFR Sgr D (see Fig. [3]). Assuming HDE316285 to be equidistant with Sgr D 
has the appeal of placing this B supergiant in the same dynamic ecosystem as the SFR 
and neighboring YSOs and maser sources, as well as a supernova remnant SNR 1.4-0.1 a 
few arcmin to the SE. Mehringer et al. (1998) have estimated the distance of Sgr D as 8 
kpc based on H2CO and CS line kinematics measured at the VLA. At this distance, the 
stellar luminosity of HDE316285 must be corrected from log(T^/T0) = 5.5 to 6.7! This 
takes into account a small reduction to Teff (:^14kK) in the most recent line-blanketed 
models, but does very little to moderate a stellar luminosity matched only by the Pistol 
Star in the GC and r] Carina. Assuming a distance of 8 kpc, then the diameter of the 
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Figure 4. Left: Spitzer/YRS spectrum of the LMC LBV R71, with principle crystalhne sihcate 
bands and PAH features indicated. Nebular dust properties have been studied by Voors et al. 
(1999) using ISO/SWS data to 25 //m. Right: Zoom on the PAH bands in the IRS spectrum. Was 
the material comprised of crystalline dust and hydrocarbons condensed in the C-rich envelope 
of an RSG phase? 

nebula is 8.6 pc, and must have a dynamical age of at least 9300 years using an upper 
constraint of 450 km/s for the time-averaged expansion velocity set by the unresolved 
IRS lines (A/AA c± 650). 

We cannot settle comfortably with this new distance to HDE3 16285, unfortunately. 
Blum & Damineli (1999) have presented K-band imaging of Sgr D, from which they 
argue that the SFR is located between 4-7.9 kpc (for doc = 8 kpc) and suggest Ak = 
1.9 mag from the nebular Br-7/radio flux ratio. No point sources are detected in the 
K-band image within the core of Sgr D, but focusing on sources within a few arcmin of 
the SFR, 18 of 34 sources have Ak < 1.5 mag (versus Ak = 3 for the GC itself). They 
devised a model with a uniform Ak = 2.2 mag screen of extinction at 4 kpc from the Sun, 
matching the actual star counts. HDE316285 has Ay = 8 mag {Ak ^ 0.8 mag), which 
is substantially lower, suggesting that the star is in the foreground. Indeed, University 
College London Echelle Spectrograph data taken in May 1995 show Nal D velocities of 
around +5 km/s, in agreement with a low distance to HDE316285 from comparison with 
the Brand & Blitz (1993) rotation curve using doc = 8kpc. 

In summary, we are confronted with conflicting evidence on the distance to Sgr D 
and HDE316285: on the one hand H2CO and CS radio line kinmatics place Sgr D at 8 
kpc, while near-IR modeling of the reddened star population close to SFR can satisfy 
observed star counts at half the distance. Admittedly the latter method involves many 
more approximations and ensuing uncertainties, but the optical line (Na D) kinematics 
and comparatively low reddening of HDE3 16285 taken at their face values do not favor 
a large distance. If the B supergiant star is in the foreground, between 2 and 4 kpc, 
then either a substantial amount of foreground material is present also in the immediate 
surrounding ISM to interact with the stellar outflow, or else the H2 lines are being formed 
in collisions between the present-day stellar wind and the outflow detected in our 24/im 
observations. 

Implications and Outlook. The fact that HDE3 16285 is surrounded by a bipolar nebula 
supports the notion suggested by Hillier et al. (1998) that this luminous, strong- winded B 
supergiant has recently passed through an LBV phase. This is more plausible at a lower 
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distance to the star with a commensurately lower Umit age of the nebula, (i.e., ^ 2300). 
We have also noted the lack of specific dust bands that are the signatures of crystalline 
material which normally condenses in the thick, slowly-expanding envelopes of RSGs, 
suggesting that the material forming the present-day nebula was rapidly ejected and 
could not crystallize under conditions of rapid cooling and low monomer densities. The 
explosive outfiows from rj Carina also lack crystalline silicates, though there are a number 
of dust bands in this N-rich, C-depleted environment which are yet to be properly identi- 
fied (Morris et al. 1999). In contrast, ISO spectroscopy of the LBVs AG Car, WRA 751, 
and R71 (Voors et al. 1999, 2000) reveal crystalline properties of the 0-rich dust and 
the presence of C-rich hydrocarbons. See Fig. [4] for the IRS spectrum of R71. These 
properties are more consistent with formation during a RSG phase, supplying appropri- 
ate temperature and density condensation conditions as well as surface chemistries in 
which N(C)/N(0) > 1, and fits with the rather low outfiow velocities and time-averaged 
mass-loss rates which yield the nebular masses estimated by Voors et al. 

We must be careful to consider certain potential attenuating effects on dust crys- 
tallinity when interpreting the mass-loss history. Specifically, silicates initially condensed 
amorphously in an eruption might be annealed under electron bombardment (Carrez 
et al. 2002a). Forsterite, an Mg-rich crystalline silicate detected in the aforementioned 
LBVs (see Fig.H]), can be formed at fiuences of ~ 10^^ e~/cm^ where (in the laboratory) 
particles are accelerated to several 10^ eV. This may occur in YSOs with magnetic fields. 
OB supergiants are not magnetically active (except in the rare Ofp?e cases) and there- 
fore the mechanism of electron bombardment would require acceleration by other means. 
Conversely, amorphization of initially crystalline dust may occur under heavy proton 
bombardment, which may explain the lack of crystalline silicates in the ISM (Carrez et 
al. 2002b, Brucato et al. 2004). Both effects of electron (re-) annealing and proton amor- 
phization require study applied to BSG/RSG environments in order to improve the use 
of dust grain properties as tracers of the evolutionary state of the underlying star during 
condensation. 
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